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' Very short quantum wires (quantum contacts) exhibit a conductance structure at the value of 

conductance close to 0.7 x 2e^ /h. The structure was discovered by Thomas et al (Phys. Rev. Lett. 
}^^' 77, 135 (1996)). Dependence of the structure on the longitudinal magnetic field was studied in 

^ , the same experiment. This dependence has clearly demonstrated that the effect is related to the 

electron spin. In the recent work (Sushkov, Phys. Rev. B 64, (2001)) it was suggested that the 
effect is caused by the development of the charge density wave within the contact. The many-body 
Ifartree-Fock calculation has confirmed the picture. The exchange electron-electron interaction is 
crucial for this calculation, so in agreement with experiment the effect is intrinsically related to the 
spin. However the dependence on the magnetic field has not been calculated yet. In the present 
paper the many-body Hartree-Fock calculation of conductance in the presence of the longitudinal 
' magnetic field is performed. 
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The quantized qogductance G = nG2, n = 1, 2,3, G2 = 2e^/ft., through a narrow quantum point contact was 
discovered in 1988El'EI. Thi^ quantization can be understood within a onc-diniensional (ID) non-interacting electron 
gas picture, see e.g. Ref.u. In the present work we are interested in a deviation-jfrom the integer quantization. 
This deviation, the so called "0.7 structure" has been found in experimental worksQcl. The structure is a shoulder- 
Sh ■ like feature or a narrow plateau ajL G w 0.76*2. More recent work demonstrates that there are some above barrier 
\ excitations related to the structuraj, and that the structure evolves down to G ~ 0.502 in longer quantum contacta]. 

■ Dependence of the structure on the longitudinal magnetic field was studied already in the pioneering worko. This 
O ' study clearly demonstrated that the effect is somehow related to the electron spin. There-JwrYe been numerous attempts 

, to explain the "0.7 structure" by spontaneous npa.Knctization of the ID quantum wireail3, or by the formation of a 

two-electtpn bound state with nonzero total spinoll3. However these scenarios contradict to the rigorous Lieb-Mattis 
T— I , theoremta that claims-that a ID many-body system with a potential interaction has zero spin in the ground state. 
' In the recent workEZi I have suggested that the effect is caused by the-development of the charge density wave 
' within the contact. This wave is a precursor for ID Wigner crystallizationlla. The total spin of the system withaut a 
magnetic field is zero in agreement with the Lieb-Mattis theoremE£l. The many-body Hartree-Fock calculationllZl has 
_ confirmed the picture and reproduced the "0.7 structure". The exchange electron-electron interaction is crucial for 
( <~>) . this calculation, in this sense the effect is intrinsically related to spin. However to probe explicitly the spin dependence 

■ one has to follow the path of the experimentErEl and to calculate the quantum contact conductance in the presence of 
, the longitudinal magnetic field. Such calculation is performed in the present workj— . 

■ The present work employs the same technique that has been developed in Rcf.EZl. This is the Hartree-Fock (HF) 
' method combined with the fictitious gauge field method. Howesjer there are some complications that are due to the 

M magnetic field. Let us first formulate the idea of the approachO. In independent particle approximation, i.e. in the 
' ' case of an ideal electron gas, a calculation of the conductance for a given transverse channel is straightforward 

O. G^'-^T, (1) 

■ h 



where T is the barrier transmission probability at Fermi energylUH. In case of interacting particles this formula is also 
Jv>( valid because before and after the potential barrier the density of electrons is high enough, and hence the interaction is 
5h negligible. However one cannot use the single particle description to calculate the transmission probability T because 
in the vicinity of the barrier the electron density is low, and hence the many-body effects are very important. To 
calculate the transmission probability T the following method is applied. Consider the liquid of electrons on a ID ring 
with a potential barrier somewhere on the ring. There is no current in the ground state of the system. Now let us 
apply a magnetic flux through the ring. This flux induces the electric current. Note that this is not a real magnetic 
fleld, this is a fictitious gauge field that generates the current without applying any voltage. On the one hand the 
current is related to the barrier transmission probability T . On the other hand, the current can be calculated using 
the HF method. This allows to find T with account of many-body effects. 

Now we have to repeat the same calculation with addition of the real longitudinal magnetic field B. This field 
does not influence an orbital dynamics, but it splits Fermi levels for electrons with spin up and spin down. So we 
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have to perform the HF calculation in the sector with nonzero total spin. According to the experimental datai'i 
the dependence of the conductance on the gate voltage evolves with the magnetic field. The field influences the 0.7 
structure, and at the value of the field i? « 5 — 10 T the conductance evolves to the double step function that one 
should expect for a simple spin split band. The effective value of the electron g factor depends on the width of , 
channel. For a wide channel it is close to the bulk limit g « 0.4, and for a narrow channel it is f/ ~ 1, see Refs.l 
There is only one channel in the quantum contact we are interested in, so the contact is narrow. On the other side 
the electrons incident on the contact are coming from the bulk and the degree of their polarization is determined by 
the bulk. Therefore for the estimate we take the intermediate value g ~ 0.7. Then the energy splitting corresponding 
to B = 5T is 

AS = g^iBB = 2 IQ-^eV = 2 IQ-^Eumt- (2) 

We use atomic units, so distances are measured in unites of Bohr radius, as — eh^ /me^, and energies are measured in 
units of Eurdt = ™fn //i^c^: where m is the effective electron mass and e is the dielectric constant. For experimental 
conditions of worksalD these values are the following: « 10~^y^m, Eunit ~ 10~^eT^. The energy scale @) corre- 
sponds to the temperature T 2 K. It is interesting to note that the "0.7 structure" depends on the temperature 
exactly on this scalecl^lZl. At this stage this is just an observation. In the present work we consider only the case of 
zero temperature. 

The Hamiltonian of the many body system that we consider is of the form 



H 



J2 \ ^^' + U{x,) - 2gMBB ■ Sij + ^ ^ V{x„x,), (3) 
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where indexes i and j numerates electrons, Si is the electron spin, Xi is the periodic coordinate on the ring of length 
L {0 < X < L), A = tt/2L is the fictitious gauge field, and B is the longitudinal magnetic field. The electron-electron 
Coulomb repulsion is of the form 

V{x,y)^^= \ (4) 

where a* « 2 is the effective width of the transverse diannel, and D{x,y) is the lengtij gf the shortest arc between the 
points X and y on the ring, for more details see Ref.EJ. To model the gate potentialclQ we use the following formula 
for the potential barrier 

In the present work we take I = 10 and d = 1. The plot of U{x) is shown in Fig.l by a solid line 
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FIG. 1. The solid line shows the external potential at Uo = 4. The dashed line shows the compensating potential Uc{x) 
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It was demonstrated in Ref.0 that to find the barrier transmission probabihty at Fermi energy with account of 
many-body effects one has to solve the many-body problem twice: without the barrier and with the barrier. The 
ratio of electric currents squared gives the transmission orpbability T = [Ju / JqY ■ This formula is valid without an 
external magnetic field. Repeating considerations of Rcf.EZi one can prove that with the magnetic field, i.e. with the 
spin splitting, the probability is given by 



where J|| is the electric current of electrons with spin up and spin down correspondingly. This formula gives the 
probability and hence due to eq. (|l]) it solves the many-body problem of conductance. 

To solve the many-body problem (^) we use the Hartree-Fock (HF) approximation. In the HF approximation the 
many body wave function is represented in the form of the Slater determinant of single particle orbitals (pi„{x). The 
index i shows the coordinate state of the orbital, and the index a = ±1/2 shows the spin state of the orbital. Each 
orbital obeys the equation 

hipia = eia^ia, (7) 

where tic is the single particle energy and h is the HF Hamiltonian 

hip,„{x) = { ^ 'rUeffix) - gtiBBa\ipi„{x) - ^ j ip*„{y)ip^„(y)V(x,y)dyLpja{x), (8) 

Ueff = u(x) + Y^ I WAv)?y{^^v)dy + Uc{x). 

The summations are performed over all filled orbitals. The compensating potential at this stage is zero, Uc — 0. We 
will use the compensating potential later to deal with some computational problems. 

For computations we use a finite grid. In the grid modification of the Hamiltonian (Q) the kinetic energy [p — A^ip 
is replaced by [2|(p(n)p ~ (p*{n + \)e^^'^ (p{n) — ip* {n)e~^^^ (p{n + 1)] /2/i^. Here h is the step of the grid and ip{n) is 
the wave function on the site n of the grid. The electric current corresponding to the grid Hamiltonian is 

= - E 4 [^U'^y^'^^A^ + 1) ~ + l)e-^^'V,o(")] • (9) 

j 

The current is conserved because of the gauge invariance of HF equations. 

In the computations we use the grid of 400 points and the total number of electrons N — N^+N^ = 158. Because of 
the computational time limitation it is very hard to increase substantially these numbers. The transmission probability 
versus the external potential Uq, see eq. (|5|), is shown in Fig. 2. It has been calculated at zero magnetic field, i. e. at 
= iV| = 79, for the ring of length L = 80. 
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FIG. 2. The transmission probability versus the external potential barrier height Uo- Length of the ring is L — 80. 



The Fig. 2 clearly indicates the "0.7 structure". It was discussed in detail in Ref.0 that the structure is due to the 
Coulomb exchange interaction on the barrier. 

To find the effect of the longitudinal magnetic field one has to repeat the above calculation with the magnetic field. 
In general terms it is straightforward, but there is a serious technical complication. In the above calculation the 
density of electrons outside the barrier is n « 158/80 ~ 2. The corresponding Fermi momentum is pp — ^lim « 3.5. 
The electron spectrum is discreet with energy splitting 5e ~ "ppn jL w 0.1. This is ten times larger than the splitting 
due to the magnetic field 5T. Thus the system is not sensitive to the magnetic field up to i? 50T. Certainly it 
is not a physical effect, it is a byproduct of the finite ring method used in the calculation. To resolve the problem in 
a straightforward way one has to increase the length of the riiig. by at least an order of magnitude. The density of 
electrons must be substantially higher than the critical densitjt£|. It means that the total number of electrons must 
be at least 1000-2000. This is computationally impossible. 

To resolve the problem we use the following combination of the means: l)Increase the length up to L = 260. 2)Make 
the electron-electron interaction (|^) dependent on the position on the ring: 

^ . (^a.|-i5)/.7i)%l-r5)/. + i) - (10) 

3) Introduce into eq. (^) an additional compensating potential Uc{x). The interaction Vd coincides with the Coulomb 
interaction ^^{x^y) in the interval [-15,15], and Yd vanishes outside of this interval. Since there is no interaction at 
large distances we can have a relatively low density of electrons at > 15. There is a prize for the translationally 
noninvariant interaction Yd- Unfortunately the interaction produces the unphysical Coulomb barrier. To compensate 
the unphysical barrier we introduce the compensating potential Uc{x) adjusted in such a way that in absence of the 
real potential barrier (C/q = 0, see eq. (||)) the density of electrons at the contact, —15 < a: < 15, is n « 1. The 
compensating potential Uc{x) is shown in Fig.l by the dashed line. The density of electrons n{x) obtained by solution 
of HF eqs. with Vd and Uc but without any external barrier {Uq = 0) is shown in Fig. 3. 
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FIG. 3. Self consistent density of electrons calculated with the modified Coulomb interaction Vd and with the compensating 
potential Uc- The external potential is zero, U{x) — 0. The length of the ring is L — 260. 

The energy splitting between single particle levels in this case is 5e k, 0.01. The spin flip of a single electron creates 
an additional spin up electron and an additional spin down hole, so it "costs" 2e « 0.02. This must be smaller than 
the magnetic splitting (^). So the minimal magnetic field we can consider is 5T. 

Hartree-Fock calculations for the long ring (L=260) with the modified Coulomb interaction ([l0| ) and with the 
corresponding compensating potential Uc are absolutely similar to that described before. In this case the effective 
potential Uefj, see eq. (|^), is not constant even without any external barrier (|^). However this effective potential is 
smooth and well below the Fermi energy. Therefore it does not influence the transmission probability. So in equation 
(|) for the tr ansmission probability the current Jq corresponds to the zero external potential, and the current Jjj 
corresponds to some given value of the external potential Uq. The transmission probability at zero magnetic field (i.e. 
N-^ = Ni — 79) versus the barrier height Uq is plotted in Fig. 4 by the solid line. 
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FIG. 4. The transmission probability versus the external potential barrier height Uo- This is the Hartree-Fock calculation 
on the ring of length L = 260 with the modified Coulomb interaction. The solid line corresponds to zero magnetic field. The 
long-dashed line corresponds to the field B = ST, and the dashed line corresponds to B — lOT. 

One can see that after smearing the average position of the structure is close to T « 0.7. However compared to Fig. 2 
the structure is much more pronounced. This happens because the modified Coulomb interaction ( |l0|) underestimates 
the Hartree screening, at the same time the exchange interaction responsible for the structure is not influenced by 
the modification. Underestimation of the Hartree screening leads also to the shift in the position of the step: it is at 
Uo ~ 3.5 in Fig. 4 and at Uq ^ 10.5 in Fig. 2 (full screening). However all these details are not that important, the "0.7 
structure" qualitatively and to large extent quantitatively is described. What is important is that the modification 
gives us the possibility to work with the magnetic field. The result of calculation for = 80 and = 78 is shown in 
Fig. 4 by the long dashed line. As it was discussed before it corresponds to the field B = 5T. The result for A^pj= 81 
and iV| ~ 77 is shown by the dashed line. It corresponds to the field B = lOT. In agreement with the dataoEl the 
field B ^ 5T shifts the position of the structure close to T w 0.5. Further increasing of the field practically does not 
infiuence the picture. 

It is very illustrative to perform also calculations in the Hartree approximation, i. e. calculations without exchange 
interaction. The results for i? = (solid line), B = 5T (long-dashed line), B ~ lOT (dashed line) arc plotted in Fig. 
5. 
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FIG. 5. The transmission probability versus the external potential height Uo calculated in the Hartree approximation. The 
solid line corresponds to zero magnetic field. The long-dashed line corresponds to the field B = 5T, and the dashed line 
corresponds to B — lOT. 
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As one should expect neither "0.7 structure" no magnetic field dependence appear in the Hartree calculation. 

In conclusion we have performed the many-body Hartree-Fock calculation of conductance in the presence of the 
longitudinal magnetic field. Results of the calculation are shown in Fig. 4 where the transmission probability (the 
conductance in units 2e^/ft.) is plotted versus the potential barrier height that models the gate potential. The solid 
line corresponds to zero magnetic field, the long-dashed line corresponds to the field B = 5T, and the dashed line 
corresponds to _B = lOT. There are some computational limitations related to work with the magnetic field. To 
overcome these limitations we had to modify the electron-electron Coulomb interaction at large distances. This is 
why the structure in conductance in Fig. 4 is more pronounced than the structure in the calculation without the 
modification shown in Fig. 2. Our calculation demonstrates that the magnetic field B ^ 5T shifts the position of the 
structure close to T « 0.5. Further increasing of the field practically does not influence tlie,picture. These results 
are in qualitative and to a large extent quantitative agreement with the expcriMcntal dataQO. The agreement gives 
further confirmation of the explanation of the "0.7 structure" suggested in RefJl3. 
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